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ABSTRACT 

We propose a set of measurements at a charm factory to determine the D^- 
D mixing parameters AM and AF and the strong phase difference S between 
doubly- Cabibbo-suppressed (DCS) and Cabibbo-favored (CP) neutral D de- 
cays into K^n^. The method can also be used to measure strong phase 
differences between other corresponding DCS and CF amplitudes. These 
phase differences are important for studies of D^-D mixing. 



The time development of decays of neutral D mesons to doubly- Cabibbo-suppressed 
(DCS) modes such as i^+vr" exhibits interesting behavior. At t = the only term in 
the amplitude is the direct DCS K~^7r~ term, but for t > a D^-D^ mixing 

contribution appears. The interference of this term with the DCS contribution involves 
the lifetime and mass differences of the neutral D mass eigenstates as well as the final- 
state strong phase difference 6 between the Cabibbo-favored (CF) K^tt" and 
DCS —>■ K~^'K~ decay amplitudes. 

A recent study of this process by the CLEG Collaboration [|^, when combined with 
direct measurements of the lifetime difference between CP-even and CP-odd neutral 
D mesons p|, ||, |[ ^, P], suggests that this strong phase difference may be large 0. 
The FOCUS Collaboration has also presented data in agreement with CLEO's results 

"'^To be submitted to Physics Letters B. 
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P]. Several attempts to estimate this phase theoretically 0, |T^, [TT|, Q involve various 
degrees of model-dependence. A recent proposal to measure the phase [|TB[ assumes that 
A/ = and A/ = 1 DCS D —>■ Kit amplitudes into 1 = 1/2 final states involve equal 
strong phases. In the present note we suggest an entirely assumption-free method for 
determining this phase experimentally. We show that related measurements are capable 
of determining D^-D mixing parameters at the percent level. 

A charm factory operating at the '?/'(3770) resonance produces D^D^ pairs in a state 
of definite charge-conjugation eigenvalue C = —. At slightly higher energies one can 
produce such pairs with C = + |H] . One may tag one of the neutral D mesons as a CP 



eigenstate through its decays into CP eigenmodes, such as ^^^^(Tr^, u, rj, rj', 0), K~^K~ , 
and TT+TT^. (These decays are important in a different context of studying the weak 
phase 7 in i? ^ D^K |jl5|, |l^].) The other neutral D meson must then have opposite 



CP if C{D°D ) = - and the same CP if C{D^D ) = + One measures its decay 
rate into K~tt^, which includes an interference between CF and DCS amplitudes. The 
measured rate thus is given in terms of the already-measured CF and DCS rates and 
the relative strong phase 6, which permits a determination of 6. (This phase plays an 
important role in a method for measuring j in B ^ D^K.) 



We now discuss details of this method, beginning with conventions and notation and 
ending with an estimate of the achievable precision. For the majority of our discussion 
we shall neglect CP violation in neutral D mixing and decays, which is expected to be 
very small in the Standard Model. Our master equations for correlated hadronic decays 
of D^D pairs, Eqs. ([lO|)-(|ll|) below, apply also to the case of CP violation. Towards 
the end of our study we will argue that including CP violation beyond the Standard 
Model has a negligible effect on the proposed measurement of the relative strong phase 
5. 

The mass eigenstates in the neutral D meson system may be defined as 

Di = p\D')+q\D') , 

D, ^ p\D')-q\^) , (1) 

where |pp + |gp = 1, with corresponding eigenvalues /ii,2 = ^1,2 — ^ri,2/2. Neglecting 
CP violation in D^-D^ mixing, we adopt the convention in which Di is the CP-odd 
and D2 the CP-even state, choosing further p = q = 1/ We define 

2 ' 2^2 ^ ^ 

M2-M1 Pa - Pi ^ 

x = , y = ^^ — ' A/i = /i2 - /ii = r(x - 2?/) . (3) 

The decay amplitudes are defined as 

{K-n+\D^) = Ae'^^ , {K-n+\^) = Ae'^^ , 6 = 6r-Sw , (4) 

where 6ji and 6w refer, respectively, to "right-sign" and "wrong-sign" strong phases in 
the K~TT'^ system. 
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With the above-mentioned convention for CP eigenstates, self-consistency requires us 



to take = cu, = —uc as in Ref. |TH]. The decay Di K^K then is forbidden. 



With this convention, in the SU(3) hmit, we would have A/ A = —VcdV*jVudV*g = 
-|-tan^ 9c — 0.05, where 9c is the Cabibbo angle. We thus take, more generally. A/ A = 
r > 0. 

We note that (§) implies {K^n'lD^) = — A*e*^^ which is equal, by a U-spin s ^ d 
substitution and a replacement of CKM elements, to [VcdV*^ /VudV*s\A* e^^^ . Thus 5 = 
in the SU(3) limit p|, 

We begin by considering the process studied by CLEO and FOCUS [^, H, which 
does not require a symmetric charm factory. A state which is identified at the time 
of production as a D^, e.g., by the decay D*^ tt^D^, evolves in time as D^(t) = 
D^U{t) + D°f4t), where 

fjf\ _ hp-if^it , p-iti2t] _ p-i/j-t ] cos(A^t/2) 1 . - 

^^^'^^ ~ 2^ ^ ~ \ tsm{Afit/2) j ■ 

Similarly D'^it) = 1^ f^(t) + D^f^(t). The time-dependence of a "wrong-sign" final 
state is governed by the amplitude 



{K-ir+lWit)) = e-'^*[(ir-7r+|D^)cos(A/it/2) + (if-7r+|D°)i sin(A/it/2)] 

c:^ {K-n+\D°)e-'^'e'^'^[r + ie'^rt{x -iy)/2] , (6) 

where the small-argument approximation is sufficient for the trigonometric functions 
since r,x,y -C 1. The ratio of wrong-sign to right-sign decays, as a function of time, is 
then 

' {K-TT+\D\t)) 



{K-n+\D<^{t)) 



Tt ■ ^ 

r -|- ~ w)^''^ 
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Tt 

where x' = xcos6 + ysm6, y' = ycos6 — xsin5. The suggestion of large 6 first arose 
because the CLEO analysis slightly favored negative y' in the term linear in t in the 
above equation, while the FOCUS result for y f^, the dominant one in terms of statistical 
error, favored y > 0. The results for y are summarized in Table 1. The evidence for 
y > now is slightly less than 2cr. CLEO reported y' = {—2.51i\[q ±0.3)% in an analysis 
allowing for CP violation and y' = {—2.3li\'^ ± 0.3)% when CP symmetry was assumed. 

If all three terms in Eq. (|^) (constant, linear, and quadratic in t) are measurable, and 
y is also measured, then it is possible in principle to determine x and 6 up to discrete 
ambiguities. However, most (but not all) estimates of x and y within the Standard 
Model are considerably less than a percent |2T|. If this is so, the last term in (|^) may 



be inaccessible even though there may exist evidence for the ry' term, and we need an 
independent determination of 6. This may be achieved through experiments at a charm 
factory, as we now show. 
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Table I: Results on y. A shorter lifetime for the CP-even neutral D corresponds to ?/ > 0. 



Reference 


Value (percent) 


Ax' 


E791 § 


0.8 ± 2.9 ± 1.0 


0.16 


FOCUS i 


3.42 ±1.39 ±0.74 


0.80 


Belle i 


-'-•'-'-3.5-2.1 


0.06 


CLEO § 


-1.1 ±2.5 ± 1.4 


1.18 


Belle § 


-'-•-'-'J-1.65 


(a) 


Average 


2.0 ± 1.2 


2.2 



(a) Not used in average since no systematic error was quoted. 



We consider a D^D pair produced in a state of definite charge-conjugation eigenvalue 
rjc = ±1. The initial states "^ncy namely, 



1 



± 



V2 



[D\P)D\-P)±D\P)D\-P)] 



evolve in time to 



71' 



-ifiit+t) 



{cos[A/i(t ± t)/2][D'>{p)D\~p) ± '^{p)D%-p)] 



±ism[Afi{t±i)/2][D%p)D\-p)±D\p)D\-P)]} , (9) 

Here t refers to the proper time of a state traveling along the +p direction, while t refers 
to one traveling along —p. 

We now consider decays of these correlated systems into various final states, searching 
in particular for interference effects depending on 6. In all cases we integrate with respect 
to proper time, since vertex separation in a symmetric e~^e~ "charm factory" is likely 
to be problematic. An early study of correlated D^D decays into specific flavor final 
states, assuming 6 = 0, was carried out by Bigi and Sanda More recently Xing |2^ 
has considered both time-dependent and time-integrated decays into correlated pairs of 
states, including some effects of a nonzero final state phase difference. However, he has 
not considered the cases (3~) and (3+) below, from which we propose to measure S. For 
completeness, we derive general expressions for time-integrated decay rates into a pair 
of final states /i and /2, from C = —1 and C = +1 D^D states, in agreement with ||23|| : 



+ 



^C=+l 



(/l,/2) 



+ 



+ 2Re 




(10) 
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where 

= {f,\D^){MD')±{f,\&){h\D') , (12) 
5^ = (MD0)(/2|D°)±(/i|D°)(/2|D°) . (13) 

The rate expressions simphfy if one of the states (say, /2) is a CP eigenstate S"^ with 
eigenvalue C = ±1: 

r--.y,.s,)H^..P^^^M2V , (14) 

1 - r 

1 (/i,^^)-|A5j (TTCy? ' 

where = (S'^|!D^), and we have used CP\D^) = —\D^). 

The expressions for products of amphtudes in Eqs. ([12D and (|13D can be easily gen- 
eralized allowing for CP violation: 

= ^A,A,{X2±Xi) , (16) 

q 

5^ = ^A,A,{l±X2Xi) , (17) 



where 



Here and everywhere we use a normalization in which r(D° — > /j) = 

Keeping terms up to order r'^,x^,y'^ in the expressions for rates, and assuming CP 
conservation we list the following results for various cases: 

C = -1 D^lf states: 

(1-) K-'n+{p)K-'n+{-p). 

1 



^A\x' + y') , (19) 

where A was defined in Eq. (4). This process serves to measure mixing effects when 
normalized by the one which follows. 

(2-) K-n+{p)K+Ti-{-p). 

^ A'^[l -2r^cos26 -^{x^ -y^)] . (20) 
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The dominant contribution is proportional to A"^, so by comparing this case with the 
previous one we can learn the combination x^+y^ describing mixing. Similar information 
could be gleaned from a fit to the time-distribution of a single tagged as described 
above, if x'^ + y'^ is sufficiently large, but this may not be the case. The small interference 
terms in ( pOf ) are probably unmeasurable. 

(3-) K-n+ip)Sd-p). 

r(-)(i^-7r+,5c) = A2A|J1 + Cre-*T(1 + Z/') 

^ ^2^1^(1 + 2Crcos5) . (21) 

Recall r ^ tan^^^^ ^ 0.05 in the SU(3) limit. By comparing rates with ^ = +1 final 
states such as K^K~ and C, = —1 final states such as Ks{p^, uj, 0) one can measure the 
ratio (1 + 2rcos(5)/(l — 2rcos5) and, given an independent measurement of r, one can 
obtain cos 6. 

(4-) K-7v+{p)i-{-p). 

Using a leptonic fiavor tag and defining Ai- = {£~X\D^), one finds 

T^~\K-n^,r) = A'AUl y")\ ■ (22) 

This process serves as a normalization for the one which follows describing the opposite- 
sign fiavor tag. 

(5-) K-'K+{p)l+{-p). 

T^-\K-7i+,t) = A'A]^y + ]^{x^ + y^)] , (23) 

where A^+ = {i~^X\D^) = A^-. By comparing this process with the previous one, we 
obtain -|- (x^ + ?/^)/2. This may be of interest for mixing parameters if r,x,y are of 
comparable size but, as mentioned, it is much more likely that x,y <^ r in which case 
this process can be used to measure r. 

(6-) S^{p)i^{-p). 

T^-\Sc,n = AlA\{l + y') . (24) 

Strictly speaking, the y"^ correction is one order higher in small parameters than we have 
been keeping, since A'g^ is already of order r. This process serves as a normalization for 
others. 

C = +1 D^lf states: 

(1+) K-7r+{p)K-7r+{-p). 

r^+\K-7r+,K-7r+) =4A^[r^ + ry' +-{x^ + y^)] . (25) 

8 

This expression gives information similar to that learned from the time-dependence (|^) 
in the decay of a single tagged neutral D. 
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(2+) K-7r+{p)K+7r-{-p). 



r^+){K-7c+,K+n-) = A^[l + 2r^cos26 + Ary -^{x^ -y^)] , (26) 

where y = y cos 6 + x sin 6. The correction terms are probably unmeasurable, so this 
process serves as a normahzation in comparison with the previous one. 

(3+) K-n+{p)S^{-p). 

T^+\K~n+,Sc) = A'A%\l-(re-''\\l-2(y + 3y') 

^ A'A\{l-2Crcos6){l-2Cy) . (27) 

There thus appears to be ample information to constrain cos 5 once r and y are well- 
enough known. 

(4+) K-7i+{p)r{-p). 

T^^\K~n^,n=A'AUl + 2ry~^ix'-y')] . (28) 

There is a small difference in comparison with the corresponding C = —1 case. 
(5+) K-7v+{p)i+{-p). 

T^+\K-7T+,i+)=A^A^,+ [r^ + 2ry' + ^{x^ + y^)] . (29) 

In contrast to the corresponding C = —1 case, this expression involves interference 
between mixing and DCS decay similar to Eq. (p5|). 

(6+) Sdp)i^{~p). 

T^-^\S^,n = AlAl{l-2Cy + 3y') . (30) 

This differs from the C = —1 case by a term of first order in y. 

An additional set of cases involves the detection of one flavor eigenstate with di- 
rection p and a different flavor eigenstate with direction —p. Examples are Ktt, Kp, 
and K*TT. Processes with opposite-flavor final states [such as K~7c~^{p)K^p^{—p)] serve 
for normalization. Processes with same-flavor eigenstates [such as K^7T^{p)K^ p^{—p)] 
from C = —1 D^D pairs then give rise to a leading term in the rate proportional to 
IrjC"**^* — rje~^^^ where i and j denote two different channels, is the ratio of DCS to 
CP amplitudes in the channel z, and 5i is the strong phase difference between right-sign 
and wrong-sign decays in that channel. Given three such measurements, and indepen- 
dent determinations of two of the three rj, one can solve for the three phase differences 
5i — 6j up to discrete ambiguities. Measurement of the third rj reduces those ambiguities 
by a factor of 2 and provides one constraint. There remains an overall ambiguity in 
the common sign of all phases. One can thereby check whether the quantities 6i differ 
from channel to channel. The comparison of strengths of K*tt and Kp Dalitz plot bands 
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in (DCS) D Kirn and (CF) D Knir decays already indicate the possibility of 
this difference [Q, since U-spin relations of flavor SU(3) appear to be violated in 
such cases. This method is unable to provide absolute values of any of the phases 6i, in 
contrast to those based on the examples (3~) and (S"*") above. 

At this point let us comment briefly on how our results may be modified in the 
presence of CP violation. In the Cabibbo-Kobayashi-Maskawa framework CP violation 
in neutral D mixing and decays, dominated by the first two generations, is very small 
and can be safely neglected. Extensions of the Standard Model could induce new sources 
of CP violation. The most likely sizable effect is a possible new CP violating phase, = 
aj:g{qA/pA), occuring in the interference between — D mixing and decay amplitudes 
into Kit or other hadronic states. Dependence on such a phase requires mixing and would 
affect, for instance, the term linear in t in Eq. (|^ and the ry' terms in Eqs. (^) and 
(PSp. However, it affects Eq. (0), from which cos 6 is obtained, only in terms quadratic in 
X, y, which we neglected. This is true also for possible CP violation in mixing, \q/p\ ^ 1. 

Other CP violating effects in Eq. (^) can come from direct CP violation in Cabibbo- 
favored K^n^ and in singly Cabibbo-suppressed 5*^ decays. This would 

introduce in the amplitudes A and As^ corrections of order \ {K^7i^\D ) / {K^ti^\D^)\ — 1 

and \{Sc\^) / {Sc\D'^) \ — 1, respectively. Such effects are expected to be very small in 
extensions of the Standard Model, where direct CP violation was shown to be negligible 
even in DCS decays [0. This can be checked directly by comparing and D branch- 



ing ratios into CP-conjugate states. Furthermore, these small effects can be shown to 
occur only at second order when combining the rate of Eq. ( PT| ) with its CP-conjugate. 

In order to estimate the total sample of events needed to perform a useful measure- 
ment of 5, we note that the rates for the processes of interest are given in terms of the 
already measured CF and DCS rates and the relative strong phase 5. Let us define an 
asymmetry 

r<-)(g,)-r'-)(g_) 

ri-)(s+) + r(-)(s_) ' ^ ' 

where r(~)(5'±) is a rate for a C = —1 D^D^ configuration to decay into a CP-eigenstate 
S± with direction —p and a flavor eigenstate such as K~Ti^ with direction +p [the case 
(3~) noted above)]. Eq. (^) implies a small asymmetry, A = 2rcos6. For a small 
asymmetry, a general result is that its error is approximately where N is 

the total number of events tagged with CP-even and CP-odd eigenstates. Thus we have 

A(cos(5) ~ — ^ . (32) 



The number of CP-tagged events decaying to K tt^ is related to the total number of 
pairs N{D'^lf) through N ^ 0.01A^(D°D°)i3(L>° ^ K-Tr+) ^ Ax 10-'^N{D°D^), 
since the branching-ratio-times-efficiency factor for tagging CP eigenstates is only about 
1.1% (the total branching ratio into CP eigenstates is larger than about 5% P^]). 
With r = 1.2tan2^ ~ 0.06 |l], g, one then has 

A(cos^)^ • (33) 

'NiD^D) 
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The cross section for e~^e 'i/'(3770) is about 10 nb at the peak while a foreseen 



integrated luminosity for a charm factory operating at this energy is about 3 fb ^ 
One can thus envision collecting 3 x 10^ DD pairs, of which half are charged and half 
are neutral. We are entitled to a factor of 2 for considering both K~7i^ and K'^Tr~ final 
states. We thus estimate that one may be able to reach an accuracy of about 0.07 in 
cos 5. 

The achievable error on (x^ + y'^Y^'^ is of the order of a percent. To see this, we 
note that the case (1~) mentioned above should yield about 15 x 10^(a;^ + y'^)[B{D'^ — > 
K~Tc^)]'^ ~ 2.2 X 10^(a;^ + ?/^) {K^ti'^){K^7i'^) events under the conditions mentioned 
above, giving rise to sensitivity to + at the level of about lO""*. Modest improve- 
ments will be possible by adding other final states. 

The parameter y' = y cos 5 — x sin 5 is obtained from the ry' terms occuring in the 
time-dependence (0) in the decay of a single tagged neutral Z), and in the rates 
and ( p9D measured by producing D^D pairs with C = +1. Since r > x,y, the accuracy 
of measuring y' is expected to be better than of measuring (x^ + y'^Y^'^. Once cos 5 is 
measured with the above calculated precision, separate measurements of x and y at a 
corresponding level of sensitivity may be achieved. 
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